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In pp it was shown that the Einstein equation can be derived as a local constitutive equation 
for an equilibrium spacetime thermodynamics. More recently, in the attempt to extend the same 
approach to the case of f(R) theories of gravity, it was found that a non-equilibrium setting 
is indeed required in order to fully describe both this theory as well as classical GR [2|. Here, 
elaborating on this point, we show that the dissipative character leading to a non-equilibrium 
spacetime thermodynamics is actually related — both in GR as well as in f(R) gravity — to 
non-local heat fluxes associated with the purely gravitational/internal degrees of freedom of the 
theory. In particular, in the case of GR we show that the internal entropy production term is 
identical to the so called tidal heating term of Hartle-Hawking. Similarly, for the case of f(R) 
gravity, we show that dissipative effects can be associated with the generalization of this term plus 
a scalar contribution whose presence is clearly justified within the scalar-tensor representation of 
the theory. Finally, we show that the allowed gravitational degrees of freedom can be fixed by 
the kinematics of the local spacetime causal structure, through the specific Equivalence Principle 
formulation. In this sense, the thermodynamical description seems to go beyond Einstein's theory 
as an intrinsic property of gravitation. 
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I. INTRODUCTION 

Over a decade ago, the connection between gravity, 
thermodynamics and quantum field theory, developed in 
70 's for black holes physics [3J IH |S], was strongly tight- 
ened up by Jacobson [T], who was able to derive the Ein- 
stein equations as equilibrium constitutive equations for 
spacetime, starting from the thermodynamical proper- 
ties of local causal horizons and the thermal nature of 
the Minkowski vacuum. 

Since then, this result has given great support to the 
idea that the well established black holes thermodynam- 
ics should be in fact extendable to some more general 
spacetime thermodynamics, where the fundamental in- 
gredients are a mix between causal horizons stationar- 
ity, quantum fields thermal behaviour and peculiar holo- 
graphic properties of gravity [Hj HI [5] EH EQ • 

However, more recently, it was realized that the ther- 
modynamical derivation of the Einstein equations needs 
a generalization to a non-equilibrium thermodynamical 
setting 0|T2]. This was firstly noticed in the attempt to 
reproduce the derivation a la Jacobson for a f(R) gravity 
theory, fn that case, making the horizon entropy propor- 
tional to a function of the Ricci scalar led to a break 
down of the local thermodynamical equilibrium. 

In order to recover the f(R) gravity field equations 
from the thermodynamical prescription, it was then nec- 
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essary to modify the equilibrium entropy balance of the 
system by considering some extra entropy production 
term. More surprisingly, the same problem was in fact 
pointed out, in the same works, even for GR, the main 
issue being again substantially related to the definition of 
entropy for the system. Nevertheless, a clear physical in- 
terpretation of the extra entropy production terms which 
come into play for the non-equilibrium thermodynamics 
derivation is still missing. 

In this work, following [2], we adopt the non- 
equilibrium description as the general proper setting for 
the above approach, fn such non-equilibrium setting, the 
Einstein equations (or their generalization for higher or- 
der gravity) can be still derived from a local thermody- 
namical condition, as a constitutive equation, but only as 
far as one separates the reversible and irreversible ther- 
modynamical contributions (with the reversible sector 
being associated with the Einstein field equations), fn 
this context, the general form for the irreversible/ viscous 
contributions and the constitutive equations for their co- 
efficients is identified by exploiting a fluid dynamics de- 
scription of the local causal horizon kinematics. We then 
argue that such terms, already interpreted in [S] as dis- 
sipative effects, should be identified as heat terms asso- 
ciated with the purely gravitational/internal degrees of 
freedom of the theory, fn particular, for the case of GR 
we find that the internal entropy production term is ex- 
actly the well known tidal heating term associated with 
dissipation of black hole horizons perturbations via grav- 
itational fluxes [13], HH [15] [16j [17]. Furthermore, by ap- 
plying the same approach to the f(R) theories of gravity, 
we actually find a generalization of the tidal heating effect 
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of GR together with an extra purely scalar dissipative 
contribution (as expected in this class of theories given 
their possible representation as scalar-tensor ones [15]). 

We start, in section [TTJ by reviewing and analyzing the 
operative framework for the spacetime thermodynamics 
approach. In particular, we further elaborate on the role 
of the equivalence principle (EP) in selecting gravita- 
tional theories in such a framework. We then proceed in 
section |III| to the derivation of the Einstein equation of 
state as in [TJ [T^] . In section TVl we describe the anal- 



ogy between a general fluid and the horizon congruence 
kinematics. This allows us to interpret (in a classical ir- 
reversible thermodynamics setting) the non-equilibrium 
features of the congruence as the analogues of the viscous 
dissipative terms for the fluid. In section [V] and |VlJ we re- 
derive the equations of motion for GR and f(R) gravity 
respectively, by starting from a general non-equilibrium 
entropy balance and looking at its reversible sector. We 
show that the irreversible contributions can in both cases 
be associated with the purely gravitational dissipative 
effects of the theory. In the Discussion we shall finally 
speculate on the possible lessons one can gather from 
our results and their hints towards an identification of 
the microscopic degrees of freedom of gravity. 



II. SPACETIME THERMODYNAMICS: 
FRAMEWORK 



THE 



The idea of a spacetime thermodynamics, as firstly pre- 
sented in [I], is built on the relation between the ther- 
mal character of quantum fields vacuum, as perceived 
by an accelerated Rindler observer, and the stationarity 
properties of the respective Killing horizon. By making 
use of the equivalence principle, the notion of Rindler 
frame can be used at a local level as an experimental 
setting for studying the local spacetime dynamics as the 
way the variation of the spacetime geometry follows from 
the energy variation of the matter fields. Here we re- 
view and extend the framework of the previous investiga- 
tions [TJ [T2] stressing its crucial points and its implicit 
assumptions. 



the Einstein and Strong formulations of the Equivalence 
Principle [19]) to introduce a local inertial frame (LIE). 
This is always allowed, provided one restricts oneself to a 
region of size ! 2 < 'K(p) , where lZ(p) gives the value of 
the smaller scale associated with the radius of curvature 
at p (which will be generically non zero). 

Within this region the metric will be approximately 
Minkowski, that is 



g a b = Vab + C(e 2 ), 



(1) 



the order of approximation being fixed by the local cur- 
vature. Finally, the local patch can be described via 
Riemann normal coordinates {x a }, such that p stays at 
x a = 0. On the introduced LIF, one can construct a local 




FIG. 1: A schematic view of the thermal Rindler wedge. The 
local causal horizon H plays the role of a diathermic barrier 
for the thermal Rindler wedge. The thermal system is con- 
stituted by the set of Rindler observer £ a moving along their 
isometry trajectories. For each of them the ensemble is de- 
scribed by p = Z~ 1 exp{— Hb/T) and perturbed by the energy 
SQ that flows across H. 



A. The Local Rindler Wedge 

The first step in the construction of a spacetime ther- 
mal system is the introduction of a local notion of hori- 
zon. In analogy with the global definition of a horizon as 
the boundary of the past of future null infinity, one can 
generally consider a local horizon at p, in a generic space- 
time (M,g a b), a s one side of the boundary of the past of 
a spacelike 2-surface patch B including p. Thereby, near 
p, the local horizon will be constituted by the congruence 
of null geodesies orthogonal to B, characterized by the 
past pointing tangent null vector £ a (see Fig. 1). 

With respect to the point p, one can then invoke Lo- 
cal Lorentz Invariance of spacetime (assumed in both 



Rindler frame (LRF) by the usual coordinate transfor- 
mations. 1 It is convenient here to choose x = %cos1i(?7k) 
and t — xsinh(?7/i), introducing an arbitrary rescaling 
factor n for the proper time, in order to have a clear la- 
bel for the Rindler wedge temperature in the following 
derivation. 

With this choice, in the neighborhood of p, the LRF 
metric will be 

ds 2 = -K 2 X 2 dv 2 + d X 2 + dy 2 + dz 2 , (2) 



Note that the construction of the LRF requires to fully use the 
equivalence principle, at least in its Einstein formulation, in order 
to identify geodesies motion 1201 . 
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corresponding to the action of a Lorentz boost, with ac- 
celeration a — 1/x- 

Therefore, within a small neighborhood of p, one can 
associate the boundary of the past of the patch B to a 
section of the approximate Killing horizon, centered in 
p. The future pointing approximate boost Killing vector 
X a , vanishing in p by definition, will be tangent to the 
null congruence comprising the causal horizon and will 
leave invariant the tangent plane to B at p. 

In terms of the approximate Killing vector x a , one can 
introduce a time 2 label v along the horizon null hypersur- 
face, defined by x a ^aV — 1. The Killing parameter v can 
be expressed in terms of the null congruence affine pa- 
rameter A. For a Killing horizon, the relation is generally 
given by 

A = -e~ KV , (3) 

so that the point p is located at infinite Killing parameter 
and at A = 0. 

As a consequence, one gets x a = (dX/dv)£ a , with 
(dX/dv) = — k\ and, in the same way, 

§=(^j0 = -KXB and a=(^ja=-KXa, (4) 

which gives some helpful relations between the Killing ex- 
pansion 9 and shear a and the respective affine geodesies 
quantities. 

B. Rindler Wedge Thermal Character 

Together with local flatness, the previous construction 
actually needs to further assume that the ground state 
of the fields living in the original spacetime is locally ap- 
proximated by the Minkowski vacuum. In this case, with 
respect to the approximate Killing vector flow, associated 
with the set of observers living in the wedge at x > \t\ 
(Rindler wedge), the vacuum state can be interpreted as 
an approximate thermal state, with a temperature 

T ~ T un y/—g 00 = — , (5) 

where T un = ha/2i: is the the Unruh temperature |21j 
for the Rindler observer. 

The expression in ^ shows that T stays constant 
throughout the Rindler wedge, because of the gravita- 
tional Doppler factor x associated with the Unruh tem- 
perature, and it is well defined on the horizon. Therefore, 
the thermal character of the Rindler state is effectively 



extended from the single Rindler observer to the whole 
wedge. 

As a further step, one can introduce an entropy for the 
system, via an entanglement argument. In the Rindler 
wedge, an accelerated observer can only access informa- 
tion on spacclikc slices bounded by the bifurcation plane. 
Thereby, since vacuum fluctuation between the inside 
and the outside of the wedge are correlated, she/he will 
perceive an entanglement entropy, which scales with the 
area of the local boundary and diverges with the den- 
sity of field states in the UV limit. However, with the 
introduction of an UV cut-off (gcncrically justified via 
the Planck scale quantum vibrations, the so called hori- 
zon zitterbewegung) one can make this entropy become 
actually proportional to the area, that is 

S = aA, (6) 

where the proportionality factor a can a priori depend 
on the nature of the quantum fields as well as be some 
complicate function of the position in spacetime |22U23j . 
Together with the temperature T, this notion of entropy 
allows to consider the local Rindler wedge with its Killing 
horizon as an analogue of a canonical ensemble (Gibbs 
state) bounded by a diathermic wall. 

C. Local Equilibrium Condition 

Since all the thermal information related to the Rindler 
wedge is somehow recorded on the horizon boundary |12j , 
one can define a notion of thermal equilibrium for the 
system in terms of the stationarity properties of the hori- 
zon. 3 

In this sense, the equilibrium state for the system is 
identified with the horizon cross-section at p, where the 
Killing expansion 9 and shear a vanish and the horizon 
is instantaneously stationary with respect to A, 

dS = 8{aA) = 0. (7) 

Given the relations in Q, it is clear that the equilib- 
rium state can be reached as long as the affine quantities 
(9, a) are not diverging in p. The local curvature of B 
at p determines the geometrical properties of the hori- 
zon null congruence and consequently the value of affine 
expansion and shear at p. Hence, in order to define an 
equilibrium surface one just need to require a suitably 
smooth curvature for B at p, without fixing a priori the 
values of affine expansion and shear. 

This is a very delicate point, as we will see that the 
properties of B at p, and the corresponding values of the 



2 Notice that there is no relation between the proper time r\ defined 
in the wedge and the Killing parameter v along the horizon. 
Nevertheless we need to keep the same scaling k for dimensional 
consistence. 



3 It can be proved formally that the stability of the Rindler metric 
(Lorentz invariance), hence the Rindler horizon stationarity, ac- 
tually implies the KMS conditions 1241 . responsible for the ther- 
mal character of the vacuum energy fluctuations as measured by 
a Rindler observer. 
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optical scalars of the associated null congruence, actually 
select the theories of gravity which may arise from the 
thermodynamical approach by fixing the gravitational 
degrees of freedom of the theory. In this sense, we will 
show that the choice of B can be directly related to the 
Equivalence Principle formulation, which plays a funda- 
mental role in the argument. 



III. REVIEW OF THE PREVIOUS RESULTS: 
THE GR CASE 



Moving away from the bifurcation surface at p, the 
spacetime will become dynamical and the presence of 
matter will distort the local Rindler causal structure, per- 
turbing the thermal state at the same time. 

As long as the departure from equilibrium is small 
and slow, it can be described in analogy with a quasi- 
static process where a suitably small amount of energy is 
thrown through the horizon. Then, one needs some con- 
dition to relate the spacetime geometry deformation to 
the variation of the fields energy content. In fact, thanks 
to the thermal properties of the system, this condition 
can actually be borrowed from classical equilibrium ther- 
modynamics. 

Indeed, for a slightly perturbed Gibbs state at tem- 
perature T, variations of entropy and internal energy are 
related by the Clausius relation, 



dS = 



SQ 
T 



(8) 



where the change in the mean energy is taken to be due 
to the fluxes into/from an unobservable region of space- 
time, which is perfectly thermalized by the horizon sys- 
tem, and hence it is assumed to consist entirely of heat. 
For a system undergoing a quasi-static process of energy 
exchange with the surroundings, relation ^ is nothing 
but a local thermodynamical equilibrium condition. 

Now, since the thermal Rindler wedge behaves like a 
Gibbs state, the first fundamental assumption in |T] con- 
sists in using relation ^ to relate the horizon entropy 
and the boost energy across the Rindler horizon, with T 
given by eq. 

As a second fundamental postulate, it is assumed a 
universal entropy density a per unit horizon area 8 A, 
such that 



dS = aSA, 



(9) 



thereby implicitly considering a constant UV cut-off for 
the fields, whereas in general one would have dS = 6 (a A) 
(i.e. a can be some spacetime function). With this choice, 
the changes of the entanglement entropy of the fields in 
the wedge can be effectively described in terms of geo- 
metrical variations of the horizon cross-section. 

Let us highlight here that the a = constant assump- 
tion, made in the GR derivation of [TJ [2], can be indeed 
recast as an explicit choice of a specific formulation of the 



Equivalence principle. As said, generally the UV cut-off 
a is fixed at the quantum gravity scale. This can be 
identified as the scale at which the gravitational action is 
of the order of the quantum of action h. For GR this is 
the standard Planck length l p = \J Gh/c 3 and hence it is 
directly related to the Newton constant. However, for a 
general scalar-tensor theory, i.e. a theory compatible just 
with the Einstein Equivalence principle (EEP) [IS], G is 
promoted to a spacetime field. As a consequence of this, 
one should expect that the cut-off will be generically po- 
sition dependent. In this sense, assuming dS — aSA is 
equivalent to assume the strong formulation of the Equiv- 
alence Principle (SEP) [TSj, hence effectively to allow for 
the only two SEP-compatible gravity theories: GR and 
Nordstrom gravity |25j . 

Given the assumption in ^ , a quantitative expression 
for the system entropy variation is obtained just by ap- 
plying the definition for the change of the horizon area 
in terms of the expansion rate of the null geodesies com- 
prising it, that is 



SA= I 1(9 dA, 

' H 



(10) 



with e indicating the 2-surface area element of the hori- 
zon cross-section. Moving away from the equilibrium sur- 
face at A = 0, along the null congruence, the infinitesimal 
evolution of is given by a linear expansion around its 
equilibrium value at p, up to the first order in A, 



9 +A d() 

9p + x dX 



0(X 2 ) 



(11) 



This first order coefficient will be determined as usual by 
the Raychaudhuri equation, 



d6 
dX 



Ral 



nanb 



(12) 



where \\<j\\ 2 stands for the squared congruence shear 



_ab _ 4 
<? Oab- 



by 



In this way, the entropy variation, up to (9(A ), is given 



dS = a 



[ edX 


- A ^ 


I H 





R ab e a £ b 



The mean energy variation of the thermal system, 
given by the boosted energy current flux of matter, is 
then described by the heat flux across the horizon as 



SQ= / T abX a dY. b , 

J H 



(14) 



where T ab is the matter stress energy tensor, while the 
volume element is given by dS 6 = edX£ b . The same 



4 Here we consider a vanishing twist, as the null congruence is 
taken hypcrsurface orthogonal. 
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quantity, with respect to the null congruence parameters, 
reads 



ed\ (-Ak) T ab £ a £ b 



(15) 



H 



At this point, asking for relation ^ to hold for all null 
vectors £ a , one can equate the O (A) integrands in (13) 
and ( 15 ). At the zcroth order in A, the value of heat flux 



at p is zero, hence one necessarily gets 9 P = 0. Then, to 
the first order, 



2tt 
ha 



T ab £ a £ b = (\\a\\ 2 + R ab £ a l b ) p . 



(16) 



Moreover, if one further requires that a p — 0, one is left 
with 



2tt 
ha 



Tab — Rab + $ 9ab , 



(17) 



where <& is an undetermined integration function. 

Eventually, by assuming the local energy conservation, 
that is V fc T a b = 0, applying the divergence operator on 
both sides of (TTtJ) , and using the contracted Bianchi iden- 
tity V b R ab = Tv a i?, one finally gets <& = — \R— A, hence 



2tt 
ha 



T„h — R n 



1 



Rg a b — A g ab , 



(18) 



where A is some arbitrary integration constant. Once the 
condition 



1 



AhG 



(19) 



is imposed, one can easily recognize the familiar Einstein 
equations. Noticeably, eq. (19 1 implies that the entropy 



density of the local Rindler horizon is the same as the 
one of a black hole. 

This result tells us that, given the entropy and energy 
conservation for the system, the local thermodynamical 
equilibrium condition is, in fact, equivalent to the Ein- 
stein equation for the local thermal spacetime. Further- 
more, the EP implies that the above construction can be 
done at any spacetime point p, and hence that equation 
(18 1 holds everywhere in spacetime. 



On the other hand, by allowing for some shear in p. 
that is for a p ^ 0, the Einstein equation is no more re- 
covered. As discussed in the previous section, the LRF 
is not sensitive to the exact value of the affine expansion 
and shear at p. Therefore, setting a v = is an unjusti- 
fied arbitrary choice. This was firstly realized in In 
fact, a non zero affine shear at p may change the way in 
which the equilibrium is approached by the system in the 
Killing frame. Given the relations in Q and the analogue 
of the expansion (111 for the affine shear, one realizes that 
the Killing shear falls off to zero at p as a ~ e~ 2KV when 



a vanishes, while only as a 



when a is non van- 



ishing. In this sense, for a non vanishing affine shear, the 
equilibrium approach can be considered slow enough for 
the system to be in a non-equilibrium regime. 



This argument was used in [5] to recast the thermody- 
namical derivation in a non-equilibrium setting, where 



dS > SQ/T. 



(20) 



In this new context, the Clausius law is replaced by the 
entropy balance law, 



dS = 6Q/T + diS, 



(21) 



and the extra shear term in < \16\ is then associated with 
the internal entropy production diS, generated by the 
system out of equilibrium. The internal entropy contri- 
bution, 0(A), has the form 



d,S 



~~h~ 



e A dX llcrl 



(22) 



H 



and, in analogy to the internal entropy production terms 
originating from the squared gradients of state variables, 
a universal property of systems with viscosity in non- 
equilibrium thermodynamics [26 , it is interpreted as an 
internal entropy production term due to some internal 
spacetime viscosity, with r/ = ha/ Ait. Furthermore, 
it is also noticed in [2] that by using (19 1, one gets 



rj = X/(16ttG) in agreement with the value obtained for 
the shear viscosity of the stretched horizon of a black 
hole in the so called membrane paradigm ESS ■ This 
result concludes the review of the argument described in 

HI 121 0. 



A. Some Remarks 

One can consider the local spacetime thermodynamics 
developed in [H El [12] as a useful experimental setting, 
where the local behaviour of gravity is tested in its dif- 
ferent formulations. Nevertheless, we believe that the 
argument still misses a clear understanding. 

First of all, it seems that the thermodynamical ap- 
proach cannot have a detailed control over the degrees of 
freedom of the resulting gravity theory. In particular, on 
the way they are effectively involved in the perturbation 
of the spacetime causal structure. Indeed, by allowing 
for some shear at p, the local equilibrium condition is 
formally split in two parts: the £ a part, which leads to 
the Einstein equation, and the d£ a part, related to the 
shear term, which remains unexplained. In this sense, at 
the d£ a level, the local equilibrium condition is broken. 

Therefore, equilibrium thermodynamics can give a 
suitable description only under the assumption that the 
affine congruence orthogonal to B has zero expansion and 
shear at p. However, we saw that this is equivalent to re- 
quire that the chosen B (and hence its associated null 
congruence) is less general than the one allowed by the 
assumed entropy-area relation (or alternatively by the 
SEP). Such an ansatz seems too restrictive to consider 
( 18 1 as a general result. 
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On the one hand, allowing for non zero affine shear at 
p give rise to further interesting clues. In some way, the 
shear contribution in ( 16 1 brings into the entropy balance 



process a new degree of freedom, which is not fixed by 
the Ricci tensor and so has nothing to do with the local 
matter energy sources. Actually, the surface shear is gen- 
erally related to the Weyl tensor and usually associated 
with the distortion on the geodesies congruence due to a 
gravitational perturbation. 

In fact, this argument opens an issue about the absence 
so far of any role for gravitational fluxes in the system 
energy perturbation mechanism. Due to their non-local 
nature, the gravitational energy fluxes cannot be taken 
into account with a proper stress-energy tensor (SET). 
However, allowing for non-local terms, as the one in \\a\\ 2 



in (16 1, seems at odds with neglecting the role of these 



non-local energy contributions. 

On the other hand, the interpretation of the internal 
entropy contribution as a by-product of some sort of vis- 
cous work on the system, given in [5J [T5], is very reason- 
able, because the term (22 1 is actually related to some 



mechanical deformation due to the presence of shear in 
the null congruence generating the horizon. 

In this sense, such a spacetime viscosity seems natu- 
rally related to the distorsive effect of a gravitational flux, 
to be intended as a local curvature perturbation which 
is independent from the Einstein equation. This suggest 
that gravitational energy fluxes can possibly play a role 
into the total entropy balance of the system without en- 
tering into the Einstein equilibrium relation. 

Starting from these remarks, we are led to reconsider 
the thermodynamical argument of this section in a fully 
non-equilibrium setting further elaborating on the ap- 
proach taken in [3J. In particular, for the motivations 
given above, we argue that the internal entropy produc- 
tion, such as (22 1, actually indicates the presence of dis- 



sipative (irreversible) processes, to be related to the con- 
formal components of the spacetime curvature, which do 
not take part into the field equations and are associated 
with purely gravitational degrees of freedom. 



IV. INTERNAL ENTROPY IN 
NON-EQUILIBRIUM THERMODYNAMICS 

In classical non-equilibrium thermodynamics, the rate 
of change of the entropy is generally written as the sum 
of two contributions: 



dS = d e S + diS, 



(23) 



where d e S is the rate of entropy exchange with the sur- 
roundings while diS comes from the process occurring 
inside the system and is a non-negative quantity, accord- 
ingly to the second law of thermodynamics. In particu- 
lar, diS is zero for reversible (quasi-static) processes and 
positive for irreversible processes. 

The Clausius relation used for the equilibrium ap- 
proach in section III, is actually equivalent to the Clau- 



sius definition of entropy for the equilibrium system, that 
is 



d e S = SQ/T and d,S = 0, 



(24) 



as, in that case, the horizon perturbation is effectively de- 
scribed as a quasi-static process occurring in continuous 
equilibrium with the surrounding. However, this defini- 
tion does not hold true any more as irreversible processes 
come into play. 

Actually, in the non-equilibrium thermodynamical set- 
ting, the Clausius definition of entropy is generalized to 
the expression 



dS = 



6Q 
T 



6N, 



(25) 



where SQ is classically referred to as compensated heat, 
that is the heat transferred between the system and its 
surroundings, while SN, the so called uncompensated 
heat , indicates the amount of entropy associated with 
the heat which is intrinsic to the system when it under- 
goes an irreversible process. 

Let us stress that the above definition is very general, 
as it does not require either an a priori specification of 
the nature of the non-equilibrium variable, or the nature 
of the process involved. It generalizes the notion of lo- 
cal equilibrium by extending the entropy balance to the 
unbalanced contributions related to the irreversible pro- 
cesses, like dissipation (see e.g. [28]). 



The generalized Clausius relation (251 is helpful in Ol- 



der to clarify the nature of the equilibrium and non- 
equilibrium contributions defining our system entropy. 
In fact, by using the definition of the non-equilibrium 
entropy given in ( 23 1 , we can write 



d e S + d{S = 



T 



SN. 



(26) 



and identify the external and internal entropy in terms of 
the compensated and uncompensated heat, respectively 

• d e S = 6Q/T, at the reversible level, 

• diS — SN, at the irreversible level. 



With this approach, the argument described in III 



quires a straightforward interpretation. Indeed, the extra 



contribution ( 22 1 introduced by the non vanishing hori- 



zon shear is an internal entropy production term allowed 
by the most general choice of the null congruence associ- 
ated with B compatible with the area-entropy relation for 
GR (that we linked to the choice of the EP formulation) . 
Therefore, it has to be seen as a by-product of the pres- 
ence of internal/purely gravitational degrees of freedom 
of the theory which can be responsible for irreversible 
dissipative processes. 

However, in order to physically identify an internal en- 
tropy contribution diS into the general expression for the 
horizon entropy given in (13), one needs a clear under- 



standing of the relation between non-equilibrium forces 
and intrinsic spacetimes properties involved. 
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Since all the thermal information of the Rindler 
wedge vacuum is recorded on the horizon boundary |12j . 
the internal spacetime variables involved in the non- 
equilibrium process should be related to the null geodesic 
congruence kinematics around p. In this sense, a possible 
way to capture non-equilibrium features of the thermal 
system is to use the analogy between the congruence bun- 
dle comprising the horizon and a classical fluid. 



A. Null Congruence as a Continuous Medium 

Given the Price and Thorne membrane approach |27l 
[30] . the local Rindler horizon can be effectively approx- 
imated by a timelike hypersurface living just inside the 
true Rindler wedge, i.e. a timelike stretched horizon. 

With a 2+1 decomposition, the timelike congruence 
comprising the stretched horizon is formally equivalent 
to a 2+1 continuous medium (fluid) living on the space- 
like two dimensional cross section of the hypersurface and 
moving with velocity v l , defined by the unit timelike vec- 
tor tangent to the hypersurface [T2"l |2"T] . 

As the stretched horizon tends to the Rindler hori- 
zon, the unit timelike vector v l approximates the Killing 
vector x" ■ Therefore, from a purely kinematical point 
of view, we can consider the velocity gradient of the 
medium, VjVj, as the equivalent of the deformation ten- 
sor field of our congruence (extrinsic curvature), in the 
Killing frame, 



B a b = h a C hb e Xc;e, 



(27) 



where the tensor h ao is the projector onto the space- 
like cross-section transverse to the congruence (trans- 
verse metric), given by 



Kb = 9ab + Xa h + Xb k a 



(28) 



while x a is the Killing vector flow, as defined in sec- 
tion II and k is an auxiliary null vector field satisfying 
k a Xa — —1 which spans, together with x a the 2-surfacc 
horizon cross section. In both cases, the tensor can be 
decomposed into trace and symmetric trace- free part, 



V« = ^(V'v^ + Vy, 



(29) 



where V.y represents the symmetric component of V^Vj, 
Sij is the identity tensor and 



V« = ^(ViVj + VjVi) 
is the deviatoric tracclcss tensor. Similarly 



(30) 



Now, for the continuous medium, the presence of Vik 
indicates a relative motion between various parts of the 
fluid and it is responsible for the irreversible viscous 
transfer of momentum into the medium and for the con- 
sequent production of internal entropy In particular, for 
the viscous medium, the mechanical dissipation will con- 
tribute to the internal entropy density dis by an amount 5 



1 f V 

J, (p 



•PS"V«), (32) 



where P™ = p vls &y + P^ ls is the viscous pressure act- 
ing on the medium, once decomposed, respectively, in its 
bulk and trace-less components [28J. 



B. Constitutive equations 



The expression for the internal entropy in ( 32 ) is a bi 



linear form, where the external stress P V1S comes coupled 
with its conjugate strain V. Typically, V is a known 
function of the internal state variable of the medium, 
while the external forces are unknown. However, in clas- 
sical irreversible thermodynamics, under a local equilib- 
rium hypothesis, that is, as far as it is possible to assume 
that local and instantaneous relations between the ther- 
mal and mechanical properties of the system are the same 
as for uniform systems at equilibrium, the deformation 
of the velocity field and the viscous pressure (mechanical 
stress) can be related by the linear constitutive equations, 



pvis 



■CV'vp, 
2?? Vtj, 



Stoke 's law (33) 
Newton's law (34) 



where £ and r\ respectively indicates the bulk and shear 
viscosity of the medium. The dissipation coefficients are 
exactly related to the time dependence of equilibrium 
fluctuations for p vls and P™, in the framework of the 
linear response theory, via the Green-Kubo relations |28j . 

Thereby, given the linear relation introduced above, 
the viscous internal entropy density can be defined in 
terms of velocity gradient components of the medium, 
that is 



dis = £ (V\) ; 



2r) 



VII 2 , 



(35) 



~ 2 

where V 



VyV«. 



Now, given the stretched horizon approximation and 
by associating the bulk term (V p v p ) to 9 and the devia- 
toric traceless tensor V.y to a a b, one can expect a dissi- 
pativc internal entropy production term for the horizon 



B a b — -^dhab + <J a b, 



(31) 



where 9 and <r a {,, are the usual scalar expansion and shear 
of the null congruence. 



Among the several dissipative processes which can contribute to 
the internal entropy production for the medium, here we just 
consider the viscous stresses. 
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congruence, of the form 



(diSy 



2rj\\a\\ 



(36) 



The above expression provides two important insights. 
First, it identifies the congruence kincmatical quantities 
which are responsible for the irreversible viscous trans- 
fer of momentum into the system, and for the conse- 
quent production of internal entropy. In our setting, the 
favorite candidate to play the role of conjugate viscous 
stresses are the spacetime curvature perturbations. In- 
deed, the evolution of the deformation tensor B ab along 
the approximated Killing flow is driven by the equation 

X c V c B ab = KB ab - B c a B cb - h a c h b d R cedfX e X f , (37) 

whose trace and trace-free part give respectively the Ray- 
chaudhuri and shear evolution (tidal force) equations 
along the horizon |29j . In vacuum, the only source term 
for both the equations is given by the electric part of the 
Weyl tensor, C ab = h a c h b e C cdefX d X f 1M M which is 
well known to be associated with gravitational perturba- 
tions. 

Secondly, the strict correspondence between equations 



( 35 1 and ( 36 1 , seems to suggest the presence of linear 



constitutive equations, analogues of ( 33 1 , ( 34 1 ) , between 
the horizon kinematical quantities 9 and a ab and some 
external viscous stress. Therefore, by extending the fluid 
analogy argument, one should be able to relate the con- 
gruence viscosity coefficients to the local curvature fluc- 
tuations via some analogue of the Green-Kubo relations. 
We leave this for future investigations. 



V. NON-EQUILIBRIUM SPACETIME 
THERMODYNAMICS: GR 

We now have a clear way to relate the non-equilibrium 
features, arising in the thermodynamical derivation of 
the Einstein equations in section |TTT| to those kinematical 
degrees of freedom of the horizon congruence which are 
turned on by the local spacetime curvature. The hori- 
zon kinematics actually defines the intrinsic spacetime 
properties involved in the irreversible processes. 

We can then reproduce the thermodynamical deriva- 
tion of the Einstein equations with a non-equilibrium ir- 
reversible thermodynamics setting, simply by starting, 
in a quite general way, from a generic spacelike 2-surface 
patch B at p, with non vanishing 9 p and a p . 

As for the previous derivation, the strong equivalence 
principle allows us to use the entropy area relation as in 
^ , that is 



dS = aSA with a = const. 



(38) 



However, in the new setting, the Clausius relation will 



be generalized to the expression (25), where SN will now 



encode all the information about both microscopic prop- 
erties and irreversible perturbations of the system. 



The new argument starts from the same definition of 
entropy, given in (13 I, 



dS = 



idX 



H 



•-X 



R ab i a i 



(39) 

Since we are now dealing with a non-equilibrium setting, 
we expect that the entropy can be expressed as a sum of 
two different contributions dS = d e S + diS. Moreover, 
for the argument given in IV A we are able to identify the 
form of the non-equilibrium, unbalanced, entropy terms. 
Therefore, we can write 



d e S = a [ ed\{6- XR ab £ a £ b ) r 



H 



(40) 
(41) 



and separate, as previously argued, the reversible and 
irreversible levels 



(40 1 = SQ/T, at the reversible level, 



(41 1 = SN, at the irreversible level. 



From the first expression above, one has 
d e S =a [ edX (0-XR ab £ a £ b ) = 



(42) 



2tt 



edXX T nb £ a £ b = ° Q 



H 



T ' 



where the heat flux is still defined by the expression in 



(15 1. Even for the non-equilibrium setting the reversible 
heat will vanish at A — 0. Thereby, at the zero order 
in A, one deduces again 9 p — 0, while, at the first order, 
the relation R ab + $g ab = (27r/Sa)T a b, is recovered for all 
null vectors £ a . Following the previous discussion this im- 
plies, together with the conservation of the matter stress- 
energy tensor, the Einstein equations if a = (AhG)^ 1 . 
On the other hand, for the irreversible level, we have, 



in accordance with ( 22 1 



SN = diS= -< 



idXXl 



(43) 



H 



This again identifies the shear contribution as an internal 
entropy term, associating it to some irreversible dissipa- 
tive process occurring in the thermal Rindler wedge. 

To get a physical interpretation of SN with respect to 
the thermal properties of the Rindler wedge, it is helpful 
to express equation ( 43 1 in terms of the Killing horizon 



parameters. In the new frame, 



SN = diS - 



H 



edv \\a\\p > 0, 



(44) 



in accordance with the second law of thermodynamics. 
By a comparison with expression ( 36 1 , one can actually 



interpret the expression in (44 1 as the standard entropy 
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production term for a fluid with shear viscosity 77, defined 

by 



27/ a 



(45) 



that is 77 = Sa/47r, in agreement with the universal re- 
lation for the viscosity to entropy density ratio found in 
the AdS/CFT context [31] • 

While the previous discussion shows that the space- 
time thermodynamics nicely fits into a non-equilibrium 
setting, we now want to take this arguments a step fur- 
ther and ask whether the expression in ( 44 ) can be effec- 



tively related to some gravitational energy flux. 

The expression for the uncompensated heat given in 



( 44 1 quantifies the energy of the system which is effec- 



tively dissipated by the viscous process, 



TSN 



aT 



(46) 



H 



Then, by substituting a = (AhG) , from the reversible 
sector of the thermodynamical approach, the quantity in 



( 46 1 reads 



TSN 



H 



edv\\a\\ 2 p , 



(47) 



which coincides with the Hartle-Hawking formula for the 
tidal heating of a classical black hole [T31 ITIirTSI [TBI IT?] . 6 
This is a striking result as it defines the internal en- 
tropy production as a purely gravitational effect. Indeed, 
it can be associated with the work done on the horizon 
by the pcrturbativc tidal field which is described by the 
electric part of the Weyl curvature tensor. 7 Further- 
more, relation ( 45 1 suggest that such a work has to be 



seen as acted upon the internal/microscopic degrees of 
freedom of the theory rather than on macroscopic quan- 
tities (in this sense ( |47| cannot be interpreted as a stan- 
dard/reversible work term). The horizon viscosity im- 
plies that such a work will be converted into internal 
heat. Hence, the presence of the internal entropy term 
can then be directly related to the process of dissipa- 
tion via gravitational/internal degrees of freedom. In 
this sense, the irreversible sector contains the informa- 
tion about the possible activation/propagation of such 
degrees of freedom of the theory. 



6 Note that both in |17| and in 1131 the Hartle-Hawking formula ex- 
pressing the relation between the horizon area variation and the 
horizon shear is utilized with a definition of the surface gravity 
e which is half of that used in I15| . 

7 Even though the magnetic part of the Weyl tensor is actually 
necessary in order to define the Weyl curvature and the equations 
governing its propagation (Bianchi identities), this part does not 
play any direct role in determining the time derivative (evolution) 
of the congruence kinematic quantities, as it is just related to 
their spatial gradients. 



VI. NON-EQUILIBRIUM SPACETIME 
THERMODYNAMICS: f(R) 

A crucial assumption in the previous derivation was 
the validity of the SEP which allowed to consider the en- 
tropy density a as a constant. One might wonder what 
are the consequences of relaxing such an assumption in 
favor of the less restrictive Einstein Equivalence Princi- 
ple (EEP). In this case, one might generically expect that 
the entropy density is promoted to a spacetime function 
(basically because the EEP implies a spacetime depen- 
dent Newton constant). However, in the definition for 
the entanglement entropy of the Rindler wedge, this im- 
plies a possibly very complicated spacetime dependence 
for the UV cut-off. Furthermore, the specific form of 
such a cut-off is not uniquely fixed by the EEP corre- 
spondingly to the fact that the latter generically allows 
for many generalized theories of gravity. 

In order to make the argument as simple as possible, 
and following the treatment of [2], we consider here the 
specific case of f(R) gravity, which is known to be equiv- 
alent to a single field scalar-tensor theory (more precisely 
a Brans-Dicke theory with oj = and a specific potential 
for the scalar field [18 ). In this case, the UV cut-off is 
known to be proportional to some function of the curva- 
ture f (R) = f'(R) (where the prime indicates the deriva- 
tive with respect to R) , playing the role of the inverse of 
the gravitational coupling. In this case, the area entropy 
relation is known to be given by 



S = a((R)e 



(48) 



where a is still a constant (albeit a priori different from 
the one considered in the previous section). 

It is easy to see that in this case the entropy variation 
along the null congruence will be 



dS 
dX 



dX 



dX 



(49) 



where, by definition e~ 1 de/dX = 9. 

Consequently, the entropy change along the horizon 
will read j2j 



dS = 



[ edX(f 



(50) 



therefore acquiring, with respect to the previous argu- 
ment, an extra contribution f coupled to the dynamics 
of the scalar function f . (Here the dot stays for differen- 
tiation with respect to A.) 

For this reason, in order to set the instantaneous sta- 
tionarity condition at p, that is dS = 0, the affine expan- 
sion is no more a good dynamical variable. In this sense, 
it is helpful to define the quantity 6 = (6 f + f ) as a sort 
of effective expansion for the congruence. Consequently, 
the equilibrium surface for the system will be fixed by 
the condition 



0, 



(51) 
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that is 9 p = — f/f, where f = t'(R)£ a R m is generally 
nonzero. In particular, this actually provides an example 
of LRF equilibrium surface, for which 9 p is always non- 
vanishing, apart from the trivial case where f is constant, 
for which the theory will be equivalent to GR. 

From section |IV A[ we could already expect that the 
presence of the non- vanishing affine expansion would pro- 
duce a non-equilibrium contribution to the system en- 
tropy. In order to get a quantitative expression for the 
entropy change in the neighborhood of p, one again can 
consider an infinitesimal deviation of the entropy from 
its equilibrium value. 



Let us then Taylor expand the integrand in ( 50 1 around 
p up to the first order in A, that is 



where 



e = e v + A 



(f 



0(A 2 ), 



-lj2 



(52) 



(53) 



One can use the Raychaudhuri equation (12) and the 
geodesic equation £ a £ b - a — 0, to obtain the 0(X) expres- 
sion for the entropy change 



dS = a ( ed\\ [(f. o6 - tR ab )t£ b - 
Jh 



(54) 



3/2f ( 



where relation (51 ) is used to substitute f f = at p 



Now, keeping the expression in (151 for the heat flux, one 



can finally reproduce the same approach used in section 

El 

At the reversible level, the generalized Clausius rela- 
tion gives 



tRab - tab + ^gab = (2ir/ha)T ab 



(55) 



where \t is an undetermined function. With the same 
argument given in section |III| one then requires the con- 
servation of the matter stress-energy tensor and use the 
contracted Bianchi identity to write the commutator of 
the covariant derivative as 2w c .[ Qb j = R ab d c v d . In this 
way, one finds 



[flU - f„;,r" = ( \i - Df 



(56) 



and thereby 



with the identification a = (4fi,G) _1 . In [2J, the same 
result was obtained starting from the entropy balance 
relation (21), assuming a = 0, and then identifying the 



extra entropy term in 9 in the second line of ( 55 1 with a 



suitable internal entropy term. There, it was also shown 
that the alternative route of keeping the 9 2 in equation 



( 56 1 is not compatible with the conservation of the matter 



energy-momentum tensor. 

Indeed, following the previous discussion, the above 
term is expected (together with a shear dependent term) 
as an unavoidable contribution related to the irreversible 



sector of the generalized Clausius relation ( 25 1 



SN = — 



[ ed\X(af) 
Jh 



(59) 



which, as explained in section IV identifies the internal 
entropy production terms of the system. 



As expected, the internal entropy in (59 1 now shows 



contribution both from scalar and tensorial degrees of 
freedom. Indeed, by using the same argument as in the 
GR case, we again have a natural interpretation for the 
expression in ( 59 ) as the dissipative function of the sys- 
tem. 

The shear squared contribution is equivalent to the one 
found for GR, with a shear viscosity coefficient which now 
takes a factor f , 



V 



hat 

4tt ' 



(60) 



as a consequence of the UV cut-off chosen for the area 
entropy relation. 

On the other hand, the internal entropy contribution 
due to the scalar degree of freedom is now given by 



diS g = - f edAA(af)~#. 



(61) 



By making use of the kinematical analogy described in 
section |IV A| and by expressing the above equation in 
the Killing frame, one is naturally led to define the bulk 
viscosity £ as 



C _ 3af 

T ~ ~2kT' 

that is C = S hat /Air, as already found in |2J. 



(62) 



A. Gravitational dissipation in scalar-tensor 
gravity 



□f- 2 / 



(57) 



Eventually, equation (57 1, together with (55) exactly 
leads, as expected, to the field equations of f[R) gravity 

tRab - f;ab + f □ f - \ f J 9ab = ^ T ab . (58) 



In order to give a physical interpretation to (61 1, one 



can use the equivalence between f(R) and scalar-tensor 
gravity, thereby interpreting f as an effective massive 
dilaton. 

The action for f(R) gravity is given by 



S = 



Ha 

47T 



d A x, 



-gf(R) 



(63) 
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By introducing an auxiliary field <p = f(i?) 8 and as- 
suming f"(R) 7^ for all R, one can take V(<p) as the 
Legendre transform of f(R) so that R — V'(ip), thereby 
rewriting the expression in ( 63 1 as 



S 



Ha 

47T 



(Pxyf^) [<pR+ V(<p)] +S„ 



The Euler-Lagrange equations, in the Jordan frame, take 
the form 

<P (Rob - \gabR^j + (g ab a - V a V h ) <p + (65) 



equivalent to field equations given in ( 58 ) . 

In this frame, by using the relation (51), one can ex- 



press the dissipated energy coupled to the bulk and shear 
viscosity in (61 1, in terms of the auxiliary scalar field tp 



TSN = - [ ed\\(ap)T 
Jh 



— z • Z I II 

~<p (fi + \\<T\\ 



Similarly to the GR case, one expect that relation 



( 66 1 expresses some purely gravitational energy loss for 



the system, this time involving both scalar and tensorial 
fluxes. The interpretation of the term related to the shear 
is straightforward as it is clearly the generalization to a 
scalar-tensor theory of the tidal heating already obtained 
for the GR case. More problematic is the interpretation 
of the bulk viscosity (purely scalar) contribution as no 
equivalent derivation as that for the tidal heating term 
in GR has be performed (to our knowledge) for scalar- 
tensor theories of gravity. 

In this direction, as a first step, one can look for the 
effective field source terms which drive the local deforma- 
tion of the null horizon congruence. Again, we consider 
the effective expansion 9 as the suitable quantity to de- 
scribe the scalar perturbations of the horizon given that 
9 = is the condition for equilibrium. 

Let us then use equation ( 53 1 as the effective expan- 



sion rate for the horizon bundle at p. Starting from this 
equation, where now f is substituted by ip, then using 
the Raychaudhuri equation and the Jordan frame field 



equation given in ( 66 1 , one gets 



2 



a nb 



(2n/ha) ^T ah t a i 



(67) 



<p-' z v a <pv b <pe a e b . 



Now, since we arc interested only the scalar contributions 
to 6 P , we can set T a b = and \\cr\\p = 0. Then, by using 



Note that the fact that in f(R) gravity the associated scalar field 
is not a generic spacetime function but rather just of R, makes 
it possible to derive a closed system of equations without having 
to assume the equations of motion of the scalar field separately. 



the equilibrium relation in the new frame 9 — — ip (p, 
we are left with 



(64) We see now that equation (68 1 identifies the quantity 



(3/2) ip~ ip as what one might define as the gravita- 
tional energy flux associated with the solely scalar field 
degrees of freedom. 

In conclusion, we can now provide a clean interpreta- 



tion of the viscous terms in ( 66 1 as those representing the 



thermal system internal energy loss due to both scalar 
and tensor gravitational energy fluxes through the hori- 
zon. In particular, by moving to the Killing frame, (66 1 
can be rewritten as 



TSN 



1 



8ttG 



e dv 



H 



-p- 1 <p 2 + ip\\af 



(69) 



where we have set a = (4HG) 1 as required by the equa- 



(66) tions of motion ( 58 1 . The striking similarity of the above 



expression with the energy loss rate due to the gravita- 
tional radiation in scalar-tensor gravity (see e.g. equation 
(10.135) of [35]) further reinforces the above suggested 
interpretation. Furthermore, we now see that also for 



the case of f(R) gravity the relations (60 1, (62 1 plus 



the above interpretation of ( 69 1 provide a correlation 



between the transmission coefficient of the gravitational 
energy through the horizon and the horizon congruence 
viscosity. 



VII. DISCUSSION 

In non-equilibrium spacetime thermodynamics the vis- 
cous dissipative effects appear to be naturally associated 
with purely gravitational energy fluxes. Their associa- 
tion to the irreversible/dissipative sector of the theory 
strongly suggests an interpretation of their nature as, 
non-local, internal heat flows associated with the inter- 
nal spacetime degrees of freedom and clarifies why in GR 
a local, background independent, description of gravita- 
tional waves is precluded. 

Noticeably, in order to recover the field equations one 
always need to effectively isolate these dissipative con- 
tributions by neatly separating the reversible and irre- 
versible sectors of the constitutive equation. The analogy 
between the horizon null congruence and a continuous 
medium allows to recognize the natural terms related to 
the irreversible sector of the entropy balance. This effec- 
tive separation of the reversible and irreversible regimes 
is further supported, at least in GR, by the fact that the 
energy contributions occurring in the equilibrium consti- 
tutive relations have a local nature, being always related 
to the gravity field sources (Ricci curvature), whereas 
the non-equilibrium terms are intrinsically non-local and 
related to those curvature components which are inde- 
pendent from the sources distribution (Weyl curvature). 
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This actually shows that the thermodynamical derivation 
of the gravitational field equations is very general as it is 
sensitive to the whole spacetime curvature. 9 

However, a different issue is the interpretation of the 
internal entropy production terms related to dissipation 
with respect to a particular spacetime solution. While 
the association between internal entropy and allowed 
form of gravitational fluxes seems quite clear (e.g. we 
showed that the energy dissipated in GR coincides ex- 
actly the Hartlc-Hawking tidal heating term), it might 
seem however puzzling that the arbitrariness in the choice 
of B allows for non-zero shear and expansion of the 
null congruence (and hence for internal entropy produc- 
tion terms) even, for example, if one imagine to have 
performed the local Rindler wedge construction in a 
Minkowski spacetime. 

In fact, the thermodynamical approach is providing us 
just with the constitutive equations of the thermal system 
associated with local Rindler wedge, not of the spacetime 
in which the latter is constructed. The arbitrariness of 
the choice of B (and hence of the thermal system proper- 
ties) implies that such equations will at most characterize 
the structure of the gravitational theory selected by the 
entropy-area relation (the EP formulation) . In this sense 
they will not be associated with physical fluxes or cur- 
vatures of the spacetime as a whole. Hence, the possible 
presence of internal entropy terms, even when the local 
Rindler wedge is constructed in flat spacetime, does not 
imply that the latter can be seen as a system in a non- 
equilibrium state. 

Of course, one might take an alternative point of view, 
and claim that the above discussion actually shows an 
intrinsic limitation of the standard construction adopted 
here, as in [TJ |2]. In addressing this issue, a possibility 
could consist in a further characterization of the local 
Rindler wedge construction. In fact, one might choose to 
construct the 2-surface B in such a way that it will be 
sensitive to the local curvature at p and reduce to a plane 
in the flat spacetime case (i.e. B would lead generically 
to a non-zero 9 p and a p but would also reduce to the stan- 
dard bifurcation surface at p for a Rindler wedge, whose 
orthogonal null congruence has 9 p — a p = 0, in the flat 
spacetime limit). For example, this could be achieved by 
constructing B as a totally geodesic 2-dimensional space- 
like sub-manifold of the spacetime passing through p [35] . 
(That is, any geodesic passing through p and there tan- 
gent to B would have to be completely contained in B.) 
Within this alternative approach, while all the formula 
would still pertain to the thermodynamical behavior of 
the local Rindler wedge at p 10 , they would now be able 



to specialize to a specific spacetime choice and hence link 
the dynamical behavior of the wedge to the actual local 
matter-curvature content of the chosen spacetime. 

We do not see at the moment a decisive argument to 
go in one sense or the other. All in all, the whole point 
of the thermodynamical approach is not to provide an 
instrument able to reconstruct the kind of spacetime one 
is living in. Rather it is aimed to put in evidence the 
thermodynamical structure and the internal degrees of 
freedom of gravitational theories. In this sense the tradi- 
tional construction, with an arbitrary B, seems sufficient. 
We plan however to further explore this issue in future 
work. 

Another important aspect of this approach highlighted 
by this work has to do with the role that the different 
Equivalence Principle formulations play in selecting the 
possible gravitational dynamics. 

The Strong Equivalence Principle implies, for a generic 
choice of B, 8 p = 0, thereby leading to the equations of 
motion of GR (with irreversible level only corresponding 
to tensorial gravitational fluxes). 11 

The Einstein Equivalence Principle, does not fix ei- 
ther 9 nor a, thereby allowing for a generalized theory 
of gravity. For the simple case of f(R), we showed how 
the scalar degree of freedom in fact produces a dissipa- 
tive contribution which is actually associated with some 
purely gravitational scalar energy flux through the sys- 
tem boundary (and which is not part of the SET of the 
scalar field). 

The presence of an unavoidable purely scalar gravita- 
tional energy flux seems to indicate that GR and scalar- 
tensor are truly separated theories. However, the very 
same thermodynamical approach, might also suggest to 
look at such gravitational theories as different regimes of 
some more general effective description of gravity (in the 
same way as a incompressible flows can be seen as a spe- 
cial regime of a general compressible fluid). This seems 
an intriguing possibility worth further investigation as it 
might lead to a unified framework which associates differ- 
ent gravitational theories with different hydrodynamical 
regimes of the analogue flow associated with the horizon 
null congruence. 

For what regards the spacetime viscosity, the quanti- 
ties r\ and £ are found to be always related to the UV 
cut-off scale of the theory through the entropy density. 
The UV cut-off scale can be identified as the scale at 
which the gravitational action is of the order of the quan- 
tum of action H, therefore where the horizon is subject to 
some quantum fluctuations, the so called horizon zitter- 



The situation in f(R) gravity is less clear due to the presence of 
the scalar field dissipative contribution in | |69[ ). It has, however, 
to be noticed that this is not a part of the stress-energy tensor of 
the scalar field for the scalar-tensor theory equivalent to f(R). 
E.g. one can talk about dissipation only with reference to the 
local Rindler wedge as spacetime as a whole as to be seen as a 



conservative system. 

We suspect that the same formulation, together with the restric- 
tion of patches B such that the null congruence shear is zero 
at any point p (taken now as a defining property of the theory), 
may imply a restriction to the conformally invariant sector of GR 
(Einstein equations plus the condition of everywhere vanishing 
Weyl), which is basically equivalent to Nordstrom gravity. 
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bewegung. This naturally suggest to interpret the space- 
time dissipative effects as a consequence of an underly- 
ing fluctuation behaviour of spacetime at the UV cut-off 
scale, as already considered in seventies by Candelas and 
Sciama |34j and nowadays explored in different terms in 
the AdS/CFT correspondence context O [35] . In this 
sense, it is interesting to note that while in principle a 
can a priori depend on the nature of the quantum fields 
and their interactions, consistency of the gravitational 
equations derived from the thermodynamical approach, 
implies a trivial relation to the Newton constant. This 
might suggest an underlying microscopic interpretation 
of gravity along the ideas of induced gravity (see e.g. [SB] 
for a related discussion). 

Finally, the separation between reversible and irre- 
versible sectors of gravity seems to indicate that dissi- 



pation can occur just at the gravitational level, without 
any contribution from matter, however phenomena like 
the Hawking radiation, with its back-reaction on the ge- 
ometry (black hole evaporation), seems to suggest that 
this separation might be a by product of the test field 
treatment so far adopted. 
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